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THEORY 0! HEAT TRANSFER IN SMOOTH AND ROUGH PIPES* 

By G. D. Mattioli 

-Ml 

BASIC CONSIDERATIONS OF FRICTION TURBULENCE 



Heat transfer' theory in the turhulent region follows 
from the theory of turbulence described in a hook hy the 
author (reference l). Prandtl (reference 2) postulates 
momentum transfer, Taylor postulates vorticity transfer, 
while Mattioli postulates both. 

Prandtl 's attack was based on the foundations of ki- 
netic theory; the momentum transfer by the molecules ex- 
plained the viscosity. But a large difference exists be- 
tween molecular and turbulent motions. The molecules are 
small and their linear momentum is sufficient to explain 
their behavior, while the elementary units which partake 
of turbulent motion are large, and both their linear and 
angular momentum must' be considered. 

In the turbulent regime, the elementary units (which 
are not separated) will be considered as soparate enti- 
ties from a dynamical point of view. The momentum (per 
unit volume) p v (p = mass density, v = velocity_vect or ) 
must be augmented by the unit angular momentum p l" u>, 
where oj = rot v, which describes the vorticity of tho 
mean motion, and I is a length which has been introduced 
on dimensional grounds. - 

Across each surface in the turbulent regime "both lin- 
ear momentum and angular momentum are transferred, and 
each exchange is dynamically equal' to a fraction of tho 
inner stresses. ' 



Let 
ty in the 
axi s , e 
du 
dr 



P e 



x coincide with the tub 
positive x-direction, r 
the eddy diffusivity, p 

is the momentum transfer per 



axis; u is tho vel^ci 
tho radius from the tunc 
the mass density; then 

unit area and time in 
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the radial direction, 
gential shear* 



-Thi-s transfer is equal to the tan- 



'rx 



P c § (kg/m 3 ) 



(1) 



Since the vorticity transfer of each element of mass does 
not vary, the corresponding stress tensor is symmetrical; 



that is, 



= T, 



is the x component of the 



■rx 'xr* 'rx 

stress, which acts on the element in a direction perpendic- 
ular to the r-axis. If T,,„ - t yt = 0, then the tensor 



is completely defined 



'rr "xx 
The equation 



ID = _ 



du 
dr 



(2) 



is the vorticity of the mean motion. In view of the defi- 
nition of e, the difference "between the vorticity trans- 
fer per unit time in thg positive and negative directions 
of r is 2tt r p e d( I u))/dr, ,and 



_d_ 
dr 



2tt r p e 



d( I 3 m) 
dr 



dr 



is the rate of increase of the vorticity "between the radius 
r and r + dr. On a unit volume "basis 



M = — 



d 

r dr 



r p e 



d( I s up 
dr 



(3) 



is the increase of vorticity per unit volume in a unit time. 
In the steady state the density of vorticity must "be con- 
served. On the surface of each volume" element AV stresses 
must exist which cause a moment about the center of the el- 
ement of - - M A 7 irrhich destroys the point-to-point vari- 
ation of 'the vorticity. Ahout the axis of rotation of the 
element which coincides with the direction of the vortices, 
conditions are symmetrical; therefore, it is natural to 



postulate that 
tensor are: 



rr 



"rx' 



etc 



the com-oonents of the stress 



°"xx = °rr = °; 



cr 



rx 



xr 



= M/2 



This tensor is anti symmetrical . The viscosity generates 
.inner forces, the correspond ing tensor "be ing; 



f rx denotes the x-c omponent of the internal pressure act- 
ing on the surface element at right angles to the r-axis. 
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du 

"xx = V rr =0; V rx =U xr =pV^ 

where v is the ' kinematic viscosity. The components of 
the total stress tensor then are: 

n 

*h = T xx + ff xx + ^>xx = °* *** * °» *xr = P (e+u) d? 

- §; t rx - P (c +U ) 4a + H ( 4) 

2 rx dr 2 

dp 

If p = constant and - — = constant, then 

dx 

1 d(r t ) dp t 

t±- :=_£.= constant (5) 

r dr dx 

for equilibrium on the element r to r + dr , of length 
dx . From these equations follow the equations of motion. 



Employing equation (4): 
1 d 



- d * . (S) 

dx 



, . du rM 
pfe +u) r — + — 
r dr L dr 2 

1 dp 

and integrating, letting — - — = - a; 

p dx 

For r = 0, the unit -shear (left side of equation) is 
equal to zero, while on the wall (r = R) , the unit shear 

(t 0 ) "becomes p a |L 

The variable v,,, = J a |j- = is callcd the fric- 

tion velocity and can he evaluated at the wall. 

A second equation rosults from tho postulate which 
fixes M. Since tho details of turbulent mixing are not 
yet known, a formal hypothesis with respect to M must 
suffice: A vortex can be generated only at the wall, for 
fluid clomonts striking the wall at other than zero tan- 
gential velocity and sticking there must experience a de- 
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formation which will cause an. eddy to "be generated. This 
vorticity generated at the walls must be destroyed in the 
interior of the fluid by the moment -M. But what is the 
law If p I s uu is constant in the r-direction, then- M 
would necessarily be zero'. One may understandably let M 
be pr op.o'rt ianal to d(p I uu)/dr. Walls are not present 
as vorticity generators' in the case of free turbulence, 
and for this case M = 0. 



Then for our case .set 



M = \fr 



d(p I 2 uO 



dr 



(8)- 



and, together with equation (3), yields 



1 _d_ 
r dr 



[ 



r pc 



d( I s u> ) 
dr 



dr 



(9) 



Further considerations based on Prandtl's similarity 
theorem (in tho neighborhood of the wall u/v* ='universal 
function of (H - r)v*;/.v) allow \|/ and I to bo treated 
as constants, where' - . ■ 



M/ = X v*; I 2 = 2(3 V ~ 



X T * 



(10) 



and x and P are universal constants of friction-gener- 
ated turbulence. X is the Karman constant = 0.406 based 
on measurements of ETikuradse . (See -reference 3.) 8 is 
harder to establish since it depends chiefly on the veloc- 
ity 'di stfibuti on -near the wall. It would appear that 
3.46 < P. < .3.5. 

THE HYDRODYHAMIC PROBLEM 



Substitute (8) into (?) and utilize (10) yielding 

( € + v) u> - 8 ^ u" = - (11) 
- v* 2 



i ' du 
where- u 1 = 



u 



ii _- 



d e u 



dr' . .» dr" 
Integration of equation ( 9 ) yields 
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€ r u" = X v*(r u' - u + K) 



(12) 



where K is the integration constant which must he fixed 
by experiment. It is the only constant still open to ques^ 
tion on theoretical grounds. 

The 'boundary conditions which must he satisfied upon 
the integration of (ll) and (12) follow from the Prandtl 
similarity theorem, and it seems natural to assume that r ,. 
the total unit shear at the wall depends on the viscosity, 
so that one may set: 



U u 1 = ^ 



= v* 



for r = R 



The velocity u at r = R must also he introduced. The 
exchange (austausch) phenomenon at the wall is damped in 
such a thin layer that one may spy that the inner interface 
of the fluid at which turbulent motion is just apparent is 
at r = H . From similarity considerations 



u = A v* 



for 



r = R 



where A is a universal dimensi onle ss quantity which de- 
pends on .P. From experiment 



and 



A = 8 .06 



A = 7.7 



for 



for 



3 = 2, 



e = 3.5 



Subject to these bound-iry conditions the integrals of equa- 
tions (ll) and (12) yield the velocity and eddy diffusiv- 
ity distributions. 



Tor large values of He = 



2Hu ,- 



v 



(u = the mean velocity) 



equations (3) and (6) may be integrated by an approximate 
procedure, the result being correct for He — » oo . 



In the laminar sublaj^er where the velocity changes 
appr e ciably with distance, while the unit shear remains 
sensibly constant, r may be replaced by R. Let R - r = 
y, the distance from the wall, then 



f 
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If u = constant, then T\ - v*y/u. Since y = F (r) 
dF v* dF _ d a F v^ d /v» dJ\_ v.% / d^F 1 dv dF\ 



and if the dots represent der i vat i ons ■ wi th ' r e spe ct to T] , 
equations (ll) and (12) become in terms of the generalized 
variables cp(= u + ) , T) (= y + ) 



<p- » (K-u) 



For large magnitudes «f Re, (K - u) — -^g— is small 

Rv * 

compared with cp and for most fluids the influence of U 
is small, thus terms involving these two quantities will 
"be neglected.. Then, the above equations ■'became 



cp + p cp = •_! and .-^ cp = - X <P and 



(14) 
and 
(15) 



which is the same form as v = constant.. That is, the in- 
fluence of a variable v is limited to the equation (13). 

Integrals of (l4).,and (15). are fixed by the conditions 
at the wall, which in turn represent the mo-mentum transfer 
mechanism at the wall. The transfer is dissipated in the 
"thin wall layer," which is to be thought of -separately 
from the main stream as the "turbulent mixing mechanism" of 
the main stream goes "over into the pattern at the wall in 
this "thin layer." The thickness of this layer is so small 
that it may be neglected and the wall process may be 
thought of as occurring in an interface". The point y = 0 
(Tj = 0) is-. then the location of this interface (i.e., the 
inner surface of the boundary layer, where the word "inner" 
means the side toward the center) and the boundary condi- 
tions will be so applied to equations (14) and (l5). 

*5n the case of an i s othermal flow the nondimensi onal di s- 
tance from the wall is T) = v^y/u . For variable tempera- 
ture u varies considerably near the wall and this must be 
taken into account by a suitable expression. 
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The mechanism of the process in the "boundary layer 
does not need to "be known in detail, "but it suffices to 
fix the velocity and the velocity gradient at the inner 
surface of the "boundary layer. The Prandtl similarity 
theorem leads to the statement that 



cp = A 



and 



cp = 1 



for 



71 = 0 



For A = 8.06 (B = 2.46) or A = 7.7 (for B 
follows from (14) and (15) that 



and if 



*ri=y=io 



is invariable 



Vdy 



) 



y=d 



2 X 



3.5) it 

(18) 



v* = A v, 



(17) 



that is, on the inner surface of the "boundary layer the 
eddy diffusivity is not zero, which follows from the con- 
cept of the separate (from the main stream) "boundary layer. 

Equations (14) and (15) are intcgrablo for B ^ 0 

and jriold a simple integral which 'is correct for He > c° , 

if the velocity di stribut ion tho interface "between the 
laminar sublayer and the turbulent stream is utilized. 
Substituting P = 0 and A = A 0 = 9.46, those integrals 
become in parametric form (see reference 1, p. 105): 



Tl 



u 



= — ( T - in T - 1 ) + A. 
T * X ' 



d y 1 f 1 , n \ c i 



(18) 



where t = cp = 
where ir = 1 . 



v dep 



and is less than 1 except at tho wall 



v* dy 

Also' - = 1 - 1 from (l4) with B = 0. 
U T 



Sma,H values of T describe the region in which the 
laminar subl-yer transforms into the turbulent core. In 
this transition layer equation (is) may bo simpli f i ed to 
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.-. : - ■ y 

. " — = - - In T + A 0 - T)- - v = -3L (19) 

v» X • 0 X * J v XT 

e 

£ 

which are asymptotic for Re >°°- , > » in the core 

for this case. In this region one may neglect v com- 
pared to e, and equations (ll) and (12) become 

e u' = - v, 2 J; e u" = X v, (u" - ^ L ^ Ji ) (20) 

the integrals of which agree with (19) for y 0. The 
hydrodynamic part has now "been solved. The equations de- 
rived agree well with the experimental measurements of 
Nikuradse (reference 3). 



THE THERMAL PROBLEM 



P 



Let 

thermal diffusivity 

unit' heat capacity at constant pressure 



y = gp _ we-ight density 

* 

# . temperature . 

The enthalpy per unit volume is *Y c^ $ (k cal/m 3 ) and 
analogously with equation (3) 



1 d 
r 3 r 



is the increase in heat content per unit time and volume 
in annular ring of width dr. If no s our ce s . exi st 



v 8£ 1 a 

Y c» — u = — — 

? 3x r dr 



r 7 Cp ( e + a) §± 



(21) 



where u(r) and e(.r) must be obtained from hydrodynamic 
considerations. 



Equation (21) may be integrated in a manner similar to 
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fch^t employed for the hydrodynamic considerations. Fbr 
instance, the rate of heat flow normal to the wall in the 
sublayer next to the wall is postulated as constant - 
which is similar to the statement that the unit shear is 
constant through the sublayer. Then 

• . ..-(€+ a) 



Br V c p . 



or, changing variables from r t o T) (equation (13)) 

f c + |i 1 (3a) 



where q is the heat rate per unit area normal to the wall 
measured positively in the direction of wall to fluid. 

Substituting the new variable, T from equation (18) 
yields 

(- - I + -) p- = S i (h - l) (23) 



a 



For gases ^ = constant, so that this equation is direct- 
ly integrable. For liquids, even though a varies but 
little with temperature,- u cannot be considered invari- 
able. It is seen that the influence of 3, (< 0.5) in equp- 

V 

tion (23) is important only for large values of T , so 
that to effect the ' formal ' int egrat i on, & will be' set equal 

to which is representative of the liquid tempera- 

ture t)° at tne wall. The integral of equation (23) be- 
comes equal to i) 0 , for t = 1. In general, 



V * <*_* Q ) = - i (in'T + In (P ro -(P P{) -l)T)) (24) 



M- c_ g 3600 v 
\irhere Pr = = — 



a 



In the turbulent nuclear flow a, as u, may be negi- 
looted with respect to e . Then equation (21) becomes 
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. -a* - - 3 

y o - ~ r u = — 

P . ox • • - or 



where u(r)" «md e(r) are "fixed by equation (20). 

An approximate solution is now sought. The velocity 
distribution falls off so rapidly in the laminar sublayer 
thqt u may be held constant^ in the turbulent core and 
equal to the mean velocity u. Then may be set 

— = - n(x) 
dx 

which is satisfied after the quieting length. Then from 
equation (25) 

2 ) 

The loft side describes the hoat flow through the wall 
(-q) as r R; therefore, 

1 

2 



q = ; E u n(x) 



and from equation (26) 

9* - _ q_ r (2?) 



G 



dr A Op S 



which is identical with equation (20). In the core 

V c„ ■ v*. u - u, 

P ( d _ « ) = i ( 28 ) 

whore the subscript l refers to any point; that is, the 
temperature and velocity profiles are similar in this re- 
gion. The temperature ^ is fixed by the requirement 
that integral equation (28) in the buffer layer between the 
laminar sublayer and the turbulent core must become equal 
t o equat i on ( 24) . 

3, tt should be noted that the approximations are also valid 
if a is not neglected by comparison with e, and that 
for Re — s*co the thickness of the semilaminar boundary 
layer becomes vani shingly small. The above transition to 
the limit toward the inner side of the semilaminar boundary 
layer is thus identified with r >R. 
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TEMPERATURE DI STRIBUT I ON AT THE WALL 

Since the numerical value of T is small in the "buf- 
fer layer, equation (24) "becomes similar to equation (18) 
when applied here: 

which, when compared with (19), may also he written as 

UzZi t. . . i . Ao + 1 (! . la Pro ) (M , 

1 o 

Equation (28) presents a similar equation for the 
turbulent core; therefore, equation (29) yields the tem- 
perature distribution in the whole tube with the excep- 
tion of a thin boundary layer whose thickness becomes zero 

for Re £»-oo. It contains the unknown temperature ($ 0 ) 

of the fluid at the wall, which must be determined by the 
.aid of thermal .conditions at the wall. Referring to the 
velocity, the velocity increase is given as A v* , which 

occurs in the boundary layer where, for the approximation 
that ^ = 0, A = A 0 = 9.46. Also, for the temperature a 

similar change must occur in the boundary layer. Dimen- 
sional considerations lead to the relation (reference 1, 
p. 306): 

6 = £ 0 - fl w = A^ f(Pr») (30) 

' c p v * 

Here $ w is the wall temperature, A# is p constant, and 

f(Pr*) is a function of Pr which could be evaluated if 

the details of the heat transfer at the walls were known. 
As it is, it is clear that the function must.be computed 
for an appropriate mean Pr* for the boundary layer. 

Even though the present thaory will not allow the de- 
termination of f(Pr*) uniquely, its range may bo estab- 
lished by moans of a discussion of two limiting conditions, 
between which the true function must lie. 
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Say € =. 0 at the inner surface' of the "boundary layer; 
then there is no mixing of the fluid in the "boundary layer 
with that of the outer region (the core) and the heat trans- 
fer must "be "by. pure conduction. Then . - 

£J = I ; e = - : 6 
dy y c p a ' y Cp a 

where 6 is the thi.ckness of. the "boundary layer. In the 
region" near the wall u/ v * is the significant length; 
therefore a likely relation is, 

8 = A. 



from which there • f oil ows: 



.- e = a* Y -^_ p r (si) 

By' comparison with (30), f(P r ) = Pr . 

In accordance with .equation (16) ' J e |y=o is no1; z ero ; 

that is, there is turbulent transfer from the inner side 
of the "boundary layer which results in the' ' int r oduct i on of 
new fluid at this interface. This fluid stream perpendic- 
ular to the interface is proportional (-reference 1, p." 20) 
to ds/dy • which, at y = 0, is 



de 



dy 



= A v* (32) 

y=o 



bo' that the quotient of the thickness A ~~ of the "bound- 

v * 

ary layer and the magnitude of equation (32) gives the 
mean time which a fluid element spends in the "boundary 
layer: 

T = A . A v* = B 3 — 2- (33) 



The transient heating of the fluid element (which is "being 
renewed from the main stream) "by conduction for a length 
of time T under temperature difference 6 at time t = 0 
is given "by the equation (reference 4): 
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(34) 



8d cP a 

A solution of = a jg- which meets the "boundary 

dt By 

conditions where \|f is the Gaussian function 

x 



V|/(x) = 



e dx 



From equation (34) the quantity of heat which flows 
to the fluid particles per unit wall-surface area is 



0, = 



y c^ a 



p r - O 



dt = 6 V c, 



7=0 




or the average heat rate i 

0, 



q = 



e v Cp v* 



T 



a. 



so that 



6 = A d 



q. 



p * 



and "by comparison with equation (30) 

f (Pr) = ,/p7 

These calculations are qualitative, but they lead to 
a likely supposition that 



f (P r ) = Pr™ where - < m < 1 



(35) 



The problem of determining the mean Pr, is yet open. 
One could set P r * = -(Pr w + Pr 0 )/2, the mean at the wall 

and on the inner surface of the boundary layer (inner al- 
ways means toward the wall center), but the logarithmic 
mean seems more basic. 
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In the vicinity of- the wall 

y 

dy 



T] = v* / — = v* y-Cl/v) 



o 



1 1 ky 

and if — = — e Uirhich is a likely supposition), then 



V o U w V w - V 0 



1/v = 



- If ^0 I- £ 



1 a 

Since p^- = — varies with l/u mostly, approximate- 
ly there nay "be written 

^ ■ (T) " < 36 ' 

?r w Pr 0 In p— 
r o 



for the value to be used in f(P r ). 

* ale 



THE EE AT TRANSFER EQUATION 



Referring to equations (29) and (30), there results 
upon summation, using the theorem expressed "by equation 
(35): 

l^fl JL + a, Pr ,». Ao + i (! ln Pro ) (37) 

The constant A$ is yet to "be determined. Experience 
teaches (reference 5) that for Pr = const Z 1 (gases) 
the ( 4> - # w )- curves and the u-curves coincide upon appro- 
priate change of scale - which result is accomplished an- 
alytically "by placing 

A$ = A 0 = 9.4 6 
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in equatioii (37). On the othe'r hand, if this constant is 
• to have universal significance, it must hold for -all val- 
ues of Pr; if such is the case, 



ohtaining the mean over the cross section 



(•-«.) -i + A 0 (P r /-l) + 1 (l - in p ro ) 

°- (38) 



1 + I* | r A 0 (P r m -l) + - (l - — - In P r ") 

u L r * X ^ p r 0 -l T oJ_ 



ffu 

which is the equation for ^ P r . Suhstituting A 0 = 9.46, 
X = 0.4 06, and 



A = 



2H dp 



u s dx 



Vu J 



yields the expression for heat transfer 
„ - q 



V c p u U - 0 W ) 



A/8 



1 + 



'A 



9.46(P ri)i m -l) +.2.46 (l - p^ 1 zi ln p r 0 ) 



(39) 



where A is gi'ven "by (reference 1, p. 133); 



= 2.005 log Re JT - 0.806 



(40) 



an equation which is in agreement with the results of 
Uikuradse (reference 3) on smooth pipes. 

It is not known how to compute Re in equation (40) 
for noni sothermal flow. The influence of viscosity is 
limited to the laminar sublayer, on the wall side of which 
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the eddy diffusivity e is, by equation ( 16) , of the order 
of magnitude of u. Since € increases greatly with y 
as one proceeds from the wall, and since- e'+ v is the com- 
bined "austausch" coefficient, it is clear that the value 
u 0 on the wall side of the boundary layer will fix the re- 
sistance to flow. Until better data are available, it will 
suffice to set 

Ee = *e 0 = — 



into equation (40) , Utilizing 

m = 0.77 (41) 

equation (39) holds well for heating and cooling. 

Actually, equation (39) contains properties which 
are not known as the result of experiment; for instance, 
Pr o and u 0 , which are fixed by the temperature of the 

sublayer on the wall side v 0 which can be fixed by know- 
ing the rate of heat transfer. For fair values of 
Prandtl'.s modulus, or for small temperature differences 

it Q ~ ij , since for sufficiently large magnitudes of Pr 

most of the temperature drop will occur in the sublayer. 
For small P r (< 0.5), or for large temperature differ- 
ences, one may employ equation (39) and the method of 
successive approximations (two will suffice). One deter- 
mines an approximate q 0 setting $ 0 = # (mean) in the 
following equation and in equation (36): 

6 = * 0 - 4 W = 9.46 q ° Pr* m (42) 

•Y Cp v* 

With those approximate properties, the calculation is re- 
peated, using equations (39) and (40)'. Finally, the mean 
temperature across the section f in area is: 




while the mixed mean is usually 'measured 

/ n Hf 



/ udf 



NACA Technical Memorandum No. 1037 



17 



Based on the mixed mean (reference 1., pp. 129-133): 



♦ A 0 CPr,"-l> ♦ |(l - 5^ in P ro ) 



0 

where u is the mean velocity across the section. Then 



y Cp u(* 1 -d w ) 



9 .46(P rj| m „-l)+2.46 (l - p 1 1 In P ro ) 



- (43) 



APPROXIMATE EXPRESSIONS FOE PRACTICAL USE 

For numerical work, the Blasius expression is useful: 
£ = 0,0396 Ee~ 0,25 (44) 

which holds well up to Ee £ 10 6 . 

When Pr is approximately unity, one may derive 

Pr m - 1 * m (Pr-lH 1 - — In P r ~ " 1 

Pr - 1 2 

For m = 0.77 and for small temperature differences 

' Nu _ 0.0396 Re -0 ' 25 



' Re p r 1 + 1.69 Ee" 0 ' 125 (P r -l) 



(45) 



which agrees well with Prandtl except that he obtains 1.74 
for the coefficient instead' of 1.69. For . m = 0.79, agree- 
ment is perfect. The problem of more accurately determining 
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m, na tur ally - i s' open., 

For large magnitudes of Py , the term ^ ln ^ T 

may he neglected. In general, one must solve the equation: 



e = 



U 



He P r 



0.0396 Re 



-0.25 



1 + 0.199 Re, 



-0.12 5 



|^9.46(P r °' 77 -l)+2.46 (l - - 1 - In P Po ) 

(46) 



Kraussold (reference 6) has proposed an empirical ex- 
pre ssion: 

— = 0.024 He -0 ' 3 P r ~ n 
Re Pr 



(47) 



where n = 0.63 for heating, and n = 0.7 for cooling. 
It is not uninteresting to compare the two formulas. (See 
the following table.) The right side of equation (46) is 
computed for" Pp.^ = Pr &> - which will hold for small temper- 
ature differences. 





Nu/^Re P r ) 


Pr 


Equation (46) 


Sciuat i on 


(47) 




Heating 




Cooling 


1 




n = 0.63 




n = 0.7 


1 


0.00596 


0 .00380 




0.0038 0 


5 


.00155 


.00138 




.00123 


10 


.000984 


.000892 




.000759 


50 


.000313 


. 000323 




. 000246 


100 


.000187 


.000209 




.000151 



It is seen that for large Pr the results from equation 
(46) lie between those for heating and cooling* 

To illustrate the method of carrying out the computa- 
tions, we will now compute the unit conductance for water 
corresponding to a mean temperature of 50° C, a mean veloc- 
ity of -100 centimeter s 'per second through a tube of 3 cen- 
timeters diameter for a' wall temperature of 100° C. 
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As a first approximation, we will choose the proper- 
ties for d =-50°; u 0 = 0.00562 cm s /sec, Pr 0 '= 3.58. 

For « = 100° C, P fw = 1,72, then P r * = 2.43 and 

Be 0 = 53,400. From equation (46), C = 0.001696, and from 

equation (42), ® v -$ 0 = 3 1°. Therefore, d 0 = 69° C. 

The propertios corresponding to the second approxima- 
tion sro: V 0 .' = * 0.00419, ?r 0 '=.2.57; Pp w '= 1»72. 

Pr* = 2.09, Ee 0 = 71,600, and from (46) 

C = — — L = 0,001724 

V Cp u U - * w ) 

For V = 988 kg/m 3 , c p = 0.998 kcal/kg° C: 

a = = 3600X0.001721X1X988X0.998 = 6120 kcal/m 2 h° 0 

w 

ROUGH TUBES 



Tubes may be thought of as completely rough if 
> 100 (references 1 nnd 7). Actually, values of 



u 

over 60 will satisfy the requirement. k is the mean 
height of the wall asperities; still "better is the equiva- 
lent height of sand roughness. In this case, the resist- 
ance equation is quadratic where 



1 R 
—P= = 2.005 log - + 1.73 



(48) 



This equation has been e s.t ahLi she d by- oxper i merit and was 
derived by the author on the assumption that the boundary- 
layer thickness (in which the wall transfer mochanism trans 
pires) is equal to k. Applying this po stula.t e- t o the ther 
mal problem in a manner similar to that obtained for a 
smooth tube yields 



6 = A, 



c p v* 



and. -A, 



V c p v* 



a 



where the mopn time of repose of the fluid particle in the 
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"boundary layer is now proportional to k/v*. On the same 
"basis as was used previously: 



n 



p * 

For a completely rough wall, the kinematic viscosity v and 
the thermal diffusivity a may ho neglected throughout the 
whole tube, *nd in place of equation (29), the following is 
obt^ ined : 

E U - « 0 ) = — - A (50) 

q 0 v* 

where A = 7.7 for 3 = 3.5 yields the "best fit (reference 
1 , p. 277) . 

From equations (49) and (50), it follows that 



:,- n 



? (* - *w> = — + A r (-i-) - A (51) 

q V* > a ' 

Averaging (not mixed mean) over the cross section, and in- 
troducing the function factor yields 

q X/8 

(52) 



y Cp u u - *„> i + /h fA r (ill) ~ a] 




There =re no d->ta as yet against which this expression can 
"be checked. Pohl 1 s (reference 8 ) measurements extend 
kv* 

only to = 12 as the formula demands, while equation 

v kv * 

(52) is good only for — — > 60. But from these measure- 
ments, it is safe to state that the term in square brackets 
in the denomination of the right side increases with in- 
kv $ 

crease as the formula demands, while it is independ- 

ent thereof in smooth tubes. (See equation (39).) 
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SUMMARY 



The heat transfer accompanying turbulent, flow in tubes 
has been treated by a new theory of wall turbulence, and a 
formula for smooth tubes has been derived (equation (39)) 

which is asymptotic at Re It agrees very well with 

the data available to date. The formula also holds for the 
flow along a flat plate if A is based on the velocity far 
away. For rough tubes, the unit conductance is shown to be 
a function of k v* /u ; the two empirical constants (A r , n) 
which appear in equation (52) cannot yet be determined bo- 
cause of lack of experimental data. 



Translation by L. M. K. Boelter, 
University 0 f California, 
Berkeley, Calif. 
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